Comparison of coannihilation effects in low-energy MSSM 



V.A. Bednyakov 

Laboratory of Nuclear Problems, Joint Institute for Nuclear Research, 
OA ■ 141980 Dubna, Russia; E-mail: bedny@nusun.jinr.ru 

o 
o 

' Abstract. The neutralino relic density is calculated in the low-energy effective Minimal 

Supersymmetric extension of Standard Model (effMSSM). The slepton-neutralino, squark- 
neutralino and neutralino/chargino-neutralino coannihilation channels are taking into account. 
The comparative study of these coannihilations is performed and demonstrated that all of them 
give the sizable contributions to the reduction of the neutralino relic density. It is shown 
■ that the predictions for direct dark matter detection rates are not strongly affected by these 

_^ ' coannihilation channels in the effMSSM. 
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1. Introduction 



A variety of data ranging from galactic rotation curves to large scale structure formation and 
the cosmic microwave background radiation imply a significant density 0.1 < Vth 2 < 0.3 [[[]] 
of so-called cold dark matter (CDM). Here Vt = p/p c and p c is the critical closure density 
of the universe, and h is the Hubble constant in units of 100 km/sec/Mpc. It is generally 
Ph| believed that most of the CDM is made of weakly-interacting massive particles (WIMPs) 
A commonly considered candidate for the WIMP is the lightest neutralino, provided it is 
the lightest supersymmetric particle (LSP) [TJ] in the Minimal Supersymmetric extension of 
Standard Model (MSSM). In most approaches the LSP is stable due to R-parity conservation 
[Q]. The neutralino, being massive, neutral and stable, often provides a sizeable contribution to 
the relic density, which is strongly model-dependent and varies by several orders of magnitude 
over the whole allowed parameter space of the MSSM. The neutralino relic density then 
can impose stringent constraints on the parameters of the MSSM, and may have important 
consequences both for studies of SUSY at colliders and in astroparticle experiments. In light 
of this and taking into account the continuing improvements in determining the abundance of 
CDM, and other components of the Universe, which have now reached an unprecedented 
precision [|J], one needs to be able to perform an accurate enough computation of the 
WIMP relic abundance, which would allow for a reliable comparison between theory and 
observation. On this way big progress in calculations of the relic density of neutralino in 
variety of supersymmetric models has been already achieved [6-41]. 

In the early universe neutralinos existed in thermal equilibrium with the cosmic thermal 
plasma. As the universe expanded and cooled, the thermal energy is no longer sufficient to 
produce neutralinos at an appreciable rate, they decouple and their number density scales 
with co-moving volume. The sparticles significantly heavier than the LSP decouple at the 
earlier time and decay into LSPs. Nevertheless there may exist some other next-to-lightest 
sparticles (NLSPs) which are not much heavier than the stable LSP. The number densities of 
the NLSPs have only slight Boltzmann suppressions with respect to the LSP number density 
when the LSP freezes out of chemical equilibrium with the thermal bath. Therefore they 
may still be present in the thermal plasma, and NLSP-LSP and NLSP-NLSP interactions hold 
LSP in thermal equilibrium resulting with significant reduction of the LSP number density. 
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These coannihilation processes can be particularly important when the LSP-LSP annihilation 
rate itself is suppressed [[15] [14], [24] ]. In the coannihilation with the LSP can be involved any 
SUSY particle, provided its mass is almost degenerate with the mass of the LSP [fl"5l, |28|l. In 
the low-energy effective MSSM (effMSSM), where one ignores restriction from unification 
assumptions and investigates the MSSM parameter space at the weak scale [ p4| , |22| , |42| , [43| ] 
there is, in principle, no preference for the next-to-lightest SUSY particle. 

The relativistic thermal averaging formalism [[T^] was extended to include coannihilation 
processes in [fZ4J, and was implemented in the DarkSusy code p5[ ] for coannihilation of 
charginos and heavier neutralinos. In was found [[24]] that such a coannihilation significantly 
decreases the relic density. The importance of the neutralino coannihilation with sferminos 

23|], stops [TJ|, R3]] and sbottoms p3] in the 



was emphasized and investigated for sleptons [L 

so-called constrained MSSM (cMSSM) [[23], ||7p or in supergravity (mSUGRA) models Q. 

In fl41[ ] the comparative study of NCC and SLC channels, exploration of relevant 
changes in the relic density and investigation of their consequences for detection of CDM 
particles were performed in the effMSSM. This paper extends our investigations [ ]4"1] ] to 
the neutralino-stop and neutralino- sbottom coannihilations and completes our consideration 
of the subject. Therefore the calculations of neutralino relic density with inclusion of the 
all relevant coannihilation channels (NCC, SLC, SQC) can be performed in the low-energy 
effMSSM. 



2. The effMSSM approach 

As free parameters in the effMSSM, we use the gaugino mass parameters M t , M 2 the 
entries to the squark and slepton mixing matrices mj~, m-., m 2 ^, m 2 ^, m| for the 1st and 

2nd generations and m^,m^,m 2 ^,m 2 ^,m 2 ^ for the 3rd generation, respectively; the 3rd 
generation trilinear soft couplings A t , A b , A T ; the mass m A of the pseudoscalar Higgs boson, 
the Higgs superpotential parameter /i, and tan (5. To reasonably reduce the parameter space 
we assumed m% = m\ = m%, m% = m% = m% , m% = m%, m% = m\ and have fixed 

U JJ Q 1 -D (^3 ii L /13 1^3 

A b = Af = [JGp. The third gaugino mass parameter M 3 defines the mass of the gluino in 
the model and is determined by means of the GUT assumption M 2 = 0.3 M 3 . The remaining 
parameters defined our effMSSM parameter space and were scanned randomly within the 
following intervals: 

-lTeV<Mi< ITeV, -2 TeV <M 2 , /i, A t < 2 TeV, 1.5 < tan/3 < 50, 
50 GeV <M A <1 TeV, 10 GeV 2 < m|, m\, m| 8 , m| 3 < 10 6 GeV 2 . 

We have included the current experimental upper limits on sparticle masses as given by the 
Particle Data Group p5|]. The limits on the rare b — > 57 decay [ ]4^ ] have also been imposed. 
The calculations of the neutralino-nucleon cross sections, and direct detection rates follow 
the description given in |Q, . The number density is governed by the Boltzmann equation 

IS 01 

(IT) 

^ + 3tfn = -<<™>(n 2 -0 (1) 

with n either being the LSP number density if there are no other coannihilating sparticles, or 
the sum over the number densities of all coannihilation partners. The index "eq" denotes the 
corresponding equilibrium value. To solve the Boltzmann equation ([!]) one needs to evaluate 
the thermally averaged neutralino annihilation cross section (erf). Without coannihilation 
processes (av) is given as the thermal average of the LSP annihilation cross-section a xx 
times relative velocity v of the annihilating LSPs (av) = (a xx v), otherwise, it is determined 
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as (cry) = (cr c fiv), where the effective thermally averaged cross-section is obtained by 
summation over coannihilating particles [ ]16| , (24[ ] 



(O-CSV) 



E 



n eq n eq " 



(2) 



The relic density is given by Q = m x n / p c , where n denotes the nowadays number 
density of the relics. For each point in the MSSM parameter space (MSSM model) we 
have evaluated the relic density of the LSP Vth 2 under the following assumptions: ignoring 
any possibility of coannihilation (IGC), taking into account only neutralino-chargino (NCC), 
slepton (SLC), or squark (SQC) coannihilations separately, as well as including all of the 
coannihilation channels (ACC). To this end in our former code p2[ ] DarkSusy procedures of 
(a e ffv) evaluation and solution of Boltzmann equation were implemented. 
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Figure 1. Effects of squark-neutralino (SQC), slepton-neutralino (SLC), and neutralino- 
chargino(neutralino) (NCC) coannihilations in effMSSM. Panels a)-d) display ratios 
^/i 2 sqc/^/i 2 igc, fi/i 2 sLc/^ft 2 iGC, r2fr 2 Ncc/^/i 2 iGC, and fi/i 2 Acc/^/i 2 SLC for the case 
when all coannihilations are included. The maximal reduction factors for all channels (NCC, 
SQC, and SLC) are of the order of 10~ 3 . Points in circles mark cosmologically interesting 
relic density 0.1 < f2/i 2 coA < 0.3. In panel a) up-going triangles correspond to stop 
coannihilations and down-going triangles correspond to sbottom coannihilations. One can 
see that stop and sbottom equally contribute. 



Coannihilations with two-body final states that can occur between neutralinos, charginos, 
sleptons, stops and sbottoms, as long as their masses are rrii < 2m x were included. The 



Feynman amplitudes for NCC, SLC and stop coannihilations were taken from DarkSusy [ p 
[[271], an d [y35|], respectively. The amplitudes for the sbottom coannihilation were obtained 
on the basis of relevant stop amplitudes from []35|]. As in [ |4T| ] the (cr^gv) and VLh 2 were 
calculated following the relevant DarkSusy routines [p5|], which were added with with codes 
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[P7[], and p5| ] in a way that guarantees the correct inclusion of SLC and SQC. We assume 
0.1 < Vth 2 < 0.3 for the cosmologically interesting region [[]]]. 



3. Coannihilation effects in the relic density 



The general view of the reduction effect on the relic density (RD) due to SQC, SLC, NCC and 
ACC are shown in Fig. |T] as ratios £7/i 2 C OA/^ 2 iGO Here £lh 2 coA is a common notation for 
^ 2 acc> ^ 2 ncc> ^ 2 sqc o r ^ 2 slo On the basis of our sampling (20000 models tested) 
the maximum RD suppression factor for NCC and SLC channels is of the same order of about 
10~ 3 . Almost the same maximal suppression is also for squark coannihilation channels. These 
results depend on the applied experimental limits on the second-lightest neutralino, chargino 
and slepton stop and sbottom masses. The current experimental limits for m?, mp,, m^±, mf, 
and m~ b are 80-90 GeV [p3|], and therefore the critical LSP mass that enables non-negligible 
NCC, SLC, and SQC contributions is also of the same order (m x > 80 GeV). From panel 
a) of the figure one can conclude that stop (up-going triangles) and sbottom (down-going 
triangles) equally contribute to reduction of RD due to coannihilations. 




LSP mass, GeV 



Figure 2. The same as in Fig. [jj but plotted together. Here £lh 2 sQc/^h 2 iQc, 
^ 2 slc/^ 2 igcn ^ 2 ncc/^ 2 igcs an d ^ 2 acc/^ 2 Slc are marked with crosses, circles, 
dots and squares, respectively. Therefore, a square filled with a cross, circle, or dot depicts a 
model that is affected only by SQC, SLC, or NCC, respectively, and any other coannihilation 
channel gives negligible contribution. Such a situation takes place for the majority of models, 
but there are some (quite few) models, given by empty squares, for which at least two 
coannihilation channels are relevant. For example, arrows in the right side of the figure 
demonstrate how reduction of RD proceeds: SLC gives no effect (f2/i 2 sLc/^ 2 iGC = 
1), SQC reduces RD with factor ri/i 2 sQc/^^ 2 iGC ~ 0.4, and finally NCC gives main 
contribution to RD suppression £1/i 2 acc/^ 2 slg ~ ^/i 2 ncc/^ 2 igc ~ 0.04 (the square 
nearly coincides with the dot). 



The circles with symbols inside depict a some kind of "constructive" reduction, when 
due to the coannihilations the relic density falls into cosmologically interesting region 0.1 < 
Qh 2 coA < 0.3. Other points present the cases when coannihilations too strongly reduce the 
relic density. One can see that NCC plays the main role in "destructive" reduction of RD, these 
channels reduce maximal number of models form cosmologically interesting region [ p4j f|TJ]. 
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Despite the fact, in "constructive" reduction of RD all coannihilation channels contribute at 
the same strength (there are almost the same number of circled points in a)-c) panels). 
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Figure 3. Ratio ilh 2 ncc/^ 2 slc, fZ/i 2 NCc/^ 2 SQC and VLh 2 SLC / 'VLh 2 SQC versus TOjjlsp _ 
m x . Open circle indicates that the f is the NLSP, star means that the light chargino x is the 
NLSP, small filled square marks the model where the second-lightest neutralino X2 is the 
NLSP. Up-going (down-going) triangles indicates that t {b) is the NLSP. 



From Fig. ^| one can see that the reduction of RD is mainly due to the only one dominant 
coannihilation channel NCC, SQC, or SLC. The other channels of coannihilation in general 
play no role or lead only to a much smaller further reduction [|4d|]. Figure [3] shows that 
for all coannihilation channels maximal RD reduction factors (less than 0.01) occur for 
mass differences m NLSP — m LSP < 20 GeV. In contrast with NCC and SLC, SQC can 
produce the same reduction effect with larger mass difference between squark and the LSP 
(ra= — m L SP ~ 150 GeV) due to the possibility of coannihilation via strong interactions. 
For NCC and SLC channels of coannihilation, relevant effects occur if the mass difference 
between the coannihilation partner and the LSP is within 15%. It was obtained that for SQC 
the relevant effects occur if the mass difference between the coannihilating squark and the 
LSP is within 50%. 

Although other coannihilation processes (including LSP annihilation with the next-to- 
NLSP (NNLSP) and next-to-NNLSP, etc), can in principal be also open from Fig. |3| one can 
conclude that dominant coannihilation channel is defined by the type of the NLSP. If next 
neutralino \2 or chargino is me NLSP, then NCC indeed dominates. Stau f (or another 
slepton) being the NLSP indeed entails a dominant SLC effect [0T]]. The SQC dominates 
when NLSP is the stop or the sbottom. 

In Fig. |] all calculated relic densities {Vth 2 i GC , f2/i 2 S Qc, fi/i 2 sLC, ^^ 2 ncc and f2/i 2 Acc) 
are depicted in the cosmologically interesting region 0.1 < Vth 2 C oA < 0.3. There is a 
quite big amount of models (mostly with m x < 250 GeV) which are completely unaffected 
by any kind of coannihilation. When at least one of coannihilation channels is relevant, 
the RD decreases and some cosmologically unviable models with Vth 2 i GC > 0.3 enter the 
cosmologically interesting range 0.1 < fi/i 2 coA < 0.3, due to NCC (squares with a dot 
inside), SLC (squares with circles inside), SQC (squares with crosses inside), or due to joint 
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Figure 4. Illustration of the shifting of effMSSM models inside and outside the cosmologically 
interesting range 0.1 < nh 2 C oA < 0.3 due to NCC, SQC and SLC. RD nh 2 i GC , ^h 2 SQC , 
ilh 2 sLC, ^ 2 ncc an d fth 2 ACC are marked with stars, crosses, circles, small dots, and squares, 
respectively. Therefore, a superposition of all symbols corresponds to a model which is totally 
untouched by coannihilation. A star-crossed circle marks a model which is untouched by SLC 
and SQC (£1/i 2 slc = ^ 2 SQC = ^ 2 igcX but shifted down due to NCC. If the corresponding 
Qh 2 Acc (which is equal to f2/i 2 Ncc) remains within this range, it still presents in the figure 
below this star-crossed circle as an empty square with a black dot inside (see short arrow). By 
analogy, an square with a circle inside gives a model which is shifted into the region due to 
SLC only (Vlh 2 acc = fWi 2 SLC)> an d if the corresponding Q,h 2 \Qc — f2/i 2 Ncc = ^ 2 SQC is 
also in the cosmologically viable range, it is located above the symbol as an crossed star with 
a dot inside (see long arrow). A quite big amount of models is shifted out of 0.1 < flh 2 < 0.3 
due to NCC (star-crossed circles). 

contribution of NCC, SQC, or/and SLC (empty squares). There are also models which enter 
the less interesting region for LSP to be CDM (fi/i 2 CO A < 0.1). The largest amount of models 
is shifted out due to NCC (star-crossed circles), and a relatively small amount of models is 
shifted out due to SLC (crossed stars with a dot inside), SQC (circles with a star and a dot 
inside), both NCC and SLC (crossed stars). There are cosmologically interesting LSPs within 
the full mass range 20 GeV < m x < 720 GeV (Fig. |]) accessible in our scan whether or not 
coannihilation channels are included. 

4. Coannihilation effects in the detection rates 

Now we consider the influence of all coannihilation channels in question (NCC, SQC and 
SLC) on prospects for direct detection of CDM neutralinos. We compare the rate predictions 
for cosmologically interesting LSPs when the RD is evaluated with or without any of 
coannihilation channel taken into account. We have seen (Fig. |j) that the RD in most 
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Figure 5. Neutralino-proton scattering cross sections for scalar (spin-independent) interaction 
(a) and axial (spin-dependent) interaction (b). Stars, crosses, circles, small dots, and squares 
correspond to 0.1 < flh 2 iQc 7 Qh 2 sQc, ^Wi 2 slc, ^^ 2 ncc, ^ 2 acc < 0.3, respectively. 
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Figure 6. Event rate for direct neutralino detection in a 73 Ge detector. Stars, crosses, circles, 
small dots, and squares correspond to 0.1 < ilh 2 iQC, ^ 2 SQC, ^ 2 slc, f^ 2 NCC, Qh 2 Acc < 
0.3, respectively. NCC decreases the maximal event rates for models with m x > 400 GeV, 
but the models with smaller LSP mass are untouched by the coannihilations. 



models with m x < 250 GeV is untouched by SQC, SLC and NCC, mostly because the 
difference m NLSP — m x is too large to yield significant effects, therefore the corresponding 
detection rates are not influenced (depicted in the figures as square filled with a star, a 
cross and a dot simultaneously). Figure |5| shows neutralino-proton scattering cross sections 
for the scalar (spin-independent) and the axial (spin-dependent) interactions. The models 
with m x < 250 GeV are hardly affected by coannihilation, and for the majority of those 
models both neutralino-proton and neutralino-neutron scattering cross sections reach values 
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a xp < 10 17 GeV 2 with the maximal cross section of order 10 15 GeV 2 . Cosmologically 
interesting models with m x > 250 GeV were influenced by coannihilations, and the maximal 
value of the neutralino-nucleon cross-section decreases from 10 -15 GeV -2 to 5 ■ 10~ 16 GeV 2 
for the models with m x > 500 GeV. In total, independently of neglection or inclusion 
of NCC, SQC and SLC the maximal scalar scattering neutralino-nucleon cross section 
reaches 10~ 16 -10 -15 GeV -2 . The spin-dependent neutralino-nucleon cross sections are 
typically higher than the spin-independent ones, and we have found the maximal values 
at 10~ 10 GeV -2 for the axial neutralino-proton and 10~ n GeV 2 for the axial neutralino- 
neutron scattering for the models which are untouched by the coannihilations. The majority 
of cosmologically interesting models yields axial neutralino-proton scattering cross sections 
in the range 5 ■ 10~ 16 GeV 2 < a XP < 2 ■ 10~ 12 GeV -2 and axial neutralino-neutron scattering 
cross sections in the range 2 ■ 10" 16 GeV~ 2 < a xn < 8 ■ 10~ 13 GeV" 2 ©. The SQC 
contribute in reduction of the cross sections, but again not significantly. Figure |^ shows the 
expected direct detection event rates calculated for a 73 Ge detector when NCC, SQC, SLC, 
and ACC are taken into account. For models with m x < 250 GeV coannihilations of any kind 
play no role. The estimations of the event rate for models with m x > 400 GeV are decreased 
mainly due to NCC Q. 




LSPmass, GeV 

Figure 7. Neutralino-proton scattering cross sections for scalar (spin-independent) interaction. 
Expectations for GENIUS detector p7[ ] and the annual-modulation region of DAMA (shaded 
region) [Q] are also given. The maximal sensitivity of GENIUS and the region of DAMA are 
located at 40 < m x < 300 GeV. 



5. Conclusion 

The neutralino relic density (RD) is calculated taking into account slepton-neutralino (SLC), 
neutralino-chargino/neutralino (NCC), and squark-neutralino (SQC) coannihilation channels 
within the low-energy effective MSSM. The maximum factors of RD decrease due to 
NCC as well as due to SQC and SLC can reach 10~ 3 , as long as the lower experimental 
limits for m f , m t ~, m~ b , and rn x ±, are of the order of 80 GeV. SQC, NCC, SLC produce 
comparable RD reduction effects in the effMSSM. For the majority of models affected by 
coannihilations it was observed that either NCC, SQC or SLC alone produces significant 
reduction of RD while the other coannihilation channels can give considerably smaller further 
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reduction. Contrary to NCC and SLC, which produce non-negligible effect only if the 
NLSP mass is smaller than 1.15m x , for SQC the relevant NLSP mass could reach 1.50m x . 
The type of the NLSP determines the dominant coannihilation channel. In the effMSSM 
all coannihilations do not imply new cosmological limits on the mass of the LSP The 
optimistic predictions for neutralino-nucleon cross sections and LSP direct detection rates 
for cosmologically interesting models are almost untouched by these coannihilations. Only 
for large m x > 400 GeV, the respectively high values are reduced, because of corresponding 
models are ruled out from the cosmological interesting region 0.1 < fi/i 2 icc < 0.3. 

From Fig. [7] one can see that the field of maximal sensitivity of the best new-generation 
CDM detectors, like GENIUS [0, as well as the annual-modulation region of DAM A 
are located at 40 < m x < 300 GeV, where coannihilation effects are almost invisible. 
Therefore despite of obvious importance of sophisticated RD calculations including complete 
set of coannihilation channels it may happen that coannihilations will play no any role at least 
for direct detection of cold dark matter. 

This work was performed in collaboration with prof. H.V. Klapdor-Kleingrothaus and 
V.Gronewold. Author thanks Yudi Santoso for making his code available, I.V. Krivosheina 
for permanent interest to the work, the Max Planck Institut fuer Kernphysik for the hospitality 
and RFBR (Grants 00-02-17587 and 02-02-04009) for support. 
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